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ABSTRACT 

We show that rigid supersymmetry theories in four dimensions can be 
extended to give supersymmetric trace (or generalized quantum) dynamics 
theories, in which the supersymmetry algebra is represented by the general- 
ized Poisson bracket of trace supercharges, constructed from fields that form 
a trace class noncommutative graded operator algebra. In particular, super- 
symmetry theories can be turned into supersymmetric matrix models this 
way. We demonstrate our results by detailed component field calculations 
for the Wess-Zumino and the supersymmetric Yang-Mills models (the latter 
with axial gauge fixing), and then show that they are also implied by a simple 
and general superspace argument. 
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1. Introduction to Trace Dynamics 



In constructing supersymmetric field theories, one usually verifies the supersymme- 
try by doing a classical Grassmann calculation, treating the bosonic fields as classical (rather 
than operator) variables and the fermionic fields as classical Grassmann (rather than oper- 
ator Grassmann) variables. Then one quantizes by replacing the classical Poisson or Dirac 
brackets by commutators/anticommutators. We shall show in this paper that for rigid su- 
persymmetry theories, a significant generalization of this standard approach is possible, in 
which the pre-quantum bosonic and fermionic fields are respectively trace class even and 
odd grade operators, such as, for example, N x N matrices whose matrix elements are re- 
spectively the even and odd grade elements of a complex Grassmann algebra. In particular, 
our results show that rigid supersymmetry theories can be extended to give supersymmet- 
ric matrix models. The requirement that the field variables be of trace class is crucial to 
our results, since in the calculations given below, cyclic permutation of operator variables 
under the trace provides the necessary commutativity, generalizing the trivial commutativ- 
ity/anticommutativity of classical field variables, for verifying both supersymmetry of the 
Lagrangian and the closure of the supersymmetry algebra. 

Our approach is based on the trace (or generalized quantum) dynamics that we have 
proposed [1] and studied with various collaborators [2]; we in fact shall use a simplified form 
of this dynamics that becomes possible when Grassmann algebras are employed to represent 
the fermion/boson distinction. Let B\ and B2 be two N x N matrices with matrix elements 
that are even grade elements of a complex Grassmann algebra, and Tr the ordinary matrix 
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trace, which obeys the cyclic property 

TrB^ = J2( B i)mn(B 2 )nm = J^W^i)™ = TrB 2 B 1 . (la) 

m,n m,n 

Similarly, let X i an d X2 be two N x N matrices with matrix elements that are odd grade 
elements of a complex Grassmann algebra, which anticommute rather than commute, so that 
the cyclic property for these takes the form 

TrXlX2 = ^2(Xl)mn(X2)nm = ~ ^(X2)mn(Xl)mn = -1**2X1 • (16) 

m,n m,n 

The cyclic/anticyclic properties of Eqs. (la, lb) are just those assumed for the trace operation 

Tr of trace dynamics, although in Refs. [1, 2] the fermionic operators were realized as 

matrices with complex matrix elements, all of which anticommute with a grading operator 

(— 1) F which formed part of the definition of Tr. Since the use of Grassmann odd fermions 

eliminates the need for the inclusion of the (— 1) F factor *, we shall continue here to use 

the notation Tr, with the understanding that fermionic matrices or operators obey Eq. (lb) 

while bosonic matrices or operators obey Eq. (la). From Eqs. (la, b), one immediately 

derives the trilinear cyclic identities 

Tr£i[£ 2 , B 3 ] =TtB 2 [B 3 , = Tr^i, B 2 \ 

Trfli{£ 2 , B 3 } =TrB 2 {B 3 , B 1 } = TrB 3 {B u B 2 } 

(lc) 

Tr£{ X i, X2} =Tr X i[x2, B] = Tr X2 [xi, B] 

Trxi{5,X2} =Tr{ X uB} X2 = ^[XuX2]B , 
which are repeatedly used below. 

* If the (— l) F construction is combined with Grassmann odd fermions one gets the 
"supertrace" str, that obeys the cyclic property strNiN 2 = strN 2 N\ for both bosonic and 
fermionic N\ >2 . We will not use the supertrace in this article. 
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The basic observation of trace dynamics is that given the trace of a polynomial P 
constructed from noncommuting matrix or operator variables, one can define a derivative 
of the c-number TrP with respect to an operator variable O by varying and then cyclically 
permuting so that in each term the factor SO stands on the right, giving the fundamental 
definition 

STrP 

STrP = Tr— SO , (2a) 
or in the condensed notation that we shall use throughout this paper, in which P = TrP, 

6P = Tr w 6 ° ■ m 

Letting L[{g r }, {q r }] be a trace Lagrangian that is a function of the bosonic or fermionic 
operators {q r } and their time derivatives, and requiring that the trace action S = f dtL be 
stationary with respect to variations of the g r 's that preserve their bosonic or fermionic type, 
one finds [1] the operator Euler-Lagrange equations 

SL d 5~L 

Sq r dt Sq r 

Defining the momentum operator p r conjugate to q r , which is of the same bosonic or fermionic 
type as q r , by 

dq r 

the trace Hamiltonian H is defined by 

H = Tr^p r g r -L . (3b) 

r 

Performing general same-type operator variations, and using the Euler-Lagrange 



equations, we find from Eq. (3b) that the trace Hamiltonian H is a trace functional 
of the operators {q r } and {p r }, 

H = H[{q r },{p r }] , (4a) 

with the operator derivatives 

-r— = -Pr, T— = , (46) 

oq r op r 

with e r = 1(— 1) according to whether q r ,p r are bosonic (fermionic) . Letting A and B be 
two trace functions of the operators {q r } and {p r }, it is convenient to define the generalized 
Poisson bracket 

{A B} - TrV e (— — - — — \ (5a) 
^ r \8q r Sp r 8q r Sp r ) 

Then using the Hamiltonian form of the equations of motion, one readily finds that for a 

general trace functional A[{g r }, {p r }], the time derivative is given by 

|A = {A,H} ; (56) 

in particular, letting A be the trace Hamiltonian H, and using the fact that the generalized 
Poisson bracket is antisymmetric in its arguments, it follows that the time derivative of H 
vanishes. 

An important property of the generalized Poisson bracket is that it satisfies [2] the 
Jacobi identity, 

{A,{B,C}} + {C,{A,B}} + {B,{C,A}} = . (6a) 

As a consequence, if Qi and Q 2 are two conserved charges, that is if 

= ^Qi = {Qi,H} , = ^Q 2 = {Q 2 ,H}, (66) 
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then their generalized Poisson bracket {Qi,Q2} also has a vanishing generalized Poisson 
bracket with H, and is conserved. This is how we will use the trace dynamics formalism to 
get representations of the Poincare supersymmetry algebra in the following sections. 

A significant feature of trace dynamics is that, as discovered by Millard [3], the 
operator [3, 4] 

C= [Qr,Pr]- Ur,Pr} (7) 

r bosons r fermions 

is conserved by the dynamics. Making the assumption (which may presuppose taking the 
N — > oo limit) that trace dynamics is ergodic, one can then analyze [4] the statistical 
mechanics of trace dynamics for the generic case in which the conserved quantities are the 
trace Hamiltonian H and the operator C. In the analysis as given in [4], the realization of 
fermions using the (— \) F construction explicitly entered the argument in two places. The 
first was in the demonstration that trace dynamics has a generalized Liouville theorem, 
which is the foundation for a statistical treatment. It is easy to see that this demonstration 
remains valid when the fermions are realized by Grassmann odd matrices without use of 
the (— \) F construction.* The second place where the (— l) F construction played a role was 
in the issue of convergence of the partition function; because Tr(— 1) F H is indefinite even 
* In the argument at the top of p. 227 of [4], a fermionic e r = — 1 was absorbed 
through the relation e r = e m e n , with e mj „ the state grading factors introduced by (— 1) F . 
When the fermions are realized through Grassmann matrices without the (— 1) F , the state 
grading factors are unity and are absent, but the interchange of the order of derivatives with 
respect to Grassmann odd matrix elements introduces an extra minus sign in the fermionic 
case, again absorbing the fermionic e r = — 1 and showing that the deviation of the Jacobian 
from unity vanishes. 
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when the operator Hamiltonian H has a positive definite bosonic part, it was necessary in 
[4] to restrict the analysis to theories in which TrH and Tr(—1) F H both generated the same 
Hamilton equations of motion, and this led to a doubling of the complexity of the statistical 
analysis. With Grassmann fermions this problem is avoided, since for the typical models 
we are studying the bosonic part of H is a positive operator, from which TV H inherits 
good positivity properties, and so the partition function can be expected to converge. The 

canonical ensemble then takes the simple form given in Eq. (48c) of [4], 

/ 9=Z- 1 exp(-TrAC , -rH) 

r (8) 
Z= d[iexp(— TrAC — rH) , 

with d/i the invariant matrix (or operator) phase space measure provided by Liouville's the- 
orem, rather than the more complicated form given in Eq. (F.l) of [4]. (As shown in [5], this 
canonical ensemble can also be derived from the corresponding microcanonical ensemble.) 
The structure of the Ward or equipartition theorems of [4] is correspondingly simplified, and 
leads as before to the conclusion that the statistical mechanics of trace dynamics is complex 
quantum field theory, with the average of the operator C playing the role of ih. As suggested 
in [4], this means that trace dynamics behaves as a pre-quantum mechanics, in which it is 
likely that the ultraviolet divergences of quantum field theory are absent. Corrections to the 
quantum field theory approximation are expected to be of order lot, with to a characteristic 
frequency of the physics in question, and so we expect the inverse of the parameter r ap- 
pearing in the canonical ensemble of Eq. (8), which has the dimension of mass, to play a 
role analogous to that of the string tension in string theories. 

2. The Wess-Zumino Model 

We begin our discussion of component field supersymmetric models with the Wess- 
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Zumino model. We follow the notational conventions of West [6] , except that we normalize 
the fermion terms in the action differently, and we always use the Majorana representation 
for the Dirac gamma matrices. Our explicit choice of 7 matrices is given in the Appendix, 
where we discuss the properties of representation covariant 7 matrix identities that take 
a particularly simple form when expressed in Majorana representation; these will play a 
significant role in our analysis. 

We start from the trace Lagrangian 



with A, B, F, G self-adjoint bosonic N x N matrices (or operators) and with x a fermionic 
4 component column vector spinor, each spin component of which is a self-adjoint fermionic 
N x N matrix (or operator). The notation % is defined by x — X T 7°> with the transpose T 
acting only on the Dirac spinor structure, so that x T is the 4 component row vector spinor 
constructed from the same N x N matrices that appear in x, and 7 is an abbreviation for 
27°. The numerical parameters A and m are respectively the coupling constant and mass. 
Equation (9a) is identical in appearance to the usual Wess-Zumino model Lagrangian, except 
that we have explicitly symmetrized the term G{A, B}; symmetrization of the other terms 
is automatic (up to total derivatives that do not contribute to the action) by virtue of the 
cyclic property of the trace. 

Taking operator variations of Eq. (9a) by using the recipe of Eqs. (2a, b), the Euler- 




m(AF + BG- XX) 



(9a) 



\[{A 2 - B 2 )F + G{A, B} - 2 X (A - il5 B) X }) 
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Lagrange equations of Eq. (2c) take the form 

d 2 A =mF + \({A, F} + {B, G} - 2xx) 
d 2 B =mG + \(-{B, F} + {A, G} + 2i X ^x) 
-fd^x =m X + A({ A X } ~ *{B, 75X}) (96) 
F =mA + X(A 2 - B 2 ) 

G=mB + X{A,B} . 
Transforming to Hamiltonian form, the canonical momenta of Eq. (3a) are 

- • T 

Px = ~Xl =*X 

Pa =d A (10a) 

p B =d B , 
and the trace Hamiltonian is given by 

H =J d 3 xTr(±ip 2 A +p% + (VA) 2 + (V5) 2 ] - ip x fi-V X 



+I(F 2 + G 2 ) - mxx + ^P x f{(A - i lb B), X }) , 



(106) 



in which F and G are understood to be the functions of A and B given by the final two lines 
of Eq. (9b), and where we have taken care to write H so that it is manifestly symmetric in 
the identical quantities p x and ix T ■ The trace three- momentum P is given by 

P = - J d 3 xTr(p A VA + p B VB + p x V X ) , (10c) 

while the conserved operator C of Eq. (7) is given by 

C = J d 3 x([A,p A ] + [B,p B ]-{ x ,p x }) , (lOd) 

with a contraction of the spinor indices in the final term of Eq. (lOd) understood. 
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Let us now perform a supersymmetry variation of the fields given by 



5 A =e X 



5B = iej 5 x 



5 X =\[F + i lb G + -fd^A + i l5 B)]e 



(11) 



5F =eYd, X 



SG = ie^^d^X 



with e a c-number Grassmann spinor (i.e., a four component spinor, the spin components of 
which are lxl Grassmann matrices). Substituting Eq. (11) into the trace Lagrangian of 
Eq. (9a), a lengthy calculation shows that when e is constant, the variation of L vanishes. 
The calculation parallels that done in the conventional c-number Lagrangian case, except 
that the trilinear cyclic identities of Eq. (lc) are used extensively in place of commutativ- 
ity/anticommutativity of the fields, and the vanishing of the terms cubic in x is most easily 
established by using the cyclic property of the trace, which implies that 



together with the cyclic identity valid for Majorana representation 7 matrices (see the Ap- 
pendix), 



^cXaXbXd = TrtcXdXaXb = TrtcXbXdX, 



a 



(12a) 




(126) 



When e is not constant, the variation of L is given by 




(13a) 



J» = - Xf [{fd u + m)(A + i lb B) + \(A 2 - B 2 + i l5 {A, B})] 



which identifies the trace supercharge Q a as 




(136) 
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where we have again taken care to express Q a symmetrically in the identical quantities p x and 
iX T ■ It is straightforward to check, using the equations of motion and the cyclic identity, that 
Tr J M is a conserved trace supercurrent, which implies that the trace supercharge is conserved. 

We are now ready to check the closure of the supersymmetry algebra under the 
generalized Poisson bracket of Eq. (5a), which for the Hamiltonian dynamics of the Wess- 
Zumino model gives 

{Q a ,Q,} = Tr[ Tr ^ 7 + ^^-^ r ^--(a-^)] . (14a) 

There are two strategies for carrying out the considerable amount of algebra involved in 
evaluating Eq. (14a). The first is to directly rearrange into the expected form, verifying 
along the way various Majorana representation 7 matrix identities that are needed; the 
second is to first Fierz transform so as to isolate a factor of the form a T T(3, and then to 
show that this yields the expected result. We shall use the first method here, and the second 
method in discussing the supersymmetric Yang-Mills model in the next section. Proceeding 
by the first method, we find that Eq. (14) rearranges, using the cyclic identities of Eq. (lc), 
into the form 

{Q a ,Qp} = aj°Pn-ujP-P , (146) 

— * 

with H and P the trace Hamiltonian and three-momentum given above. The 7 matrix 
identities needed can be obtained by repeated applications either of the cyclic identity of 
Eq. (12b), or of the additional identity (with £, m, n spatial indices, and e^ mn the three index 
antisymmetric tensor with 6123 = 1) 

iLfcd + Idbfca ~ h e i5)ab(fl5)cd ~ tflb) dbiflb) ca 

(15) 

= 5ad(fldbc - {fll)adhc + £tmn{iamlb)ad{ian)cb , 
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which we have verified by the method described in the Appendix. It is also easy to check 
that Q e plays the role of the generator of supersymmetry transformations for the dynamical 
variables A, B, x under the generalized Poisson bracket, since we readily find (for constant 
Grassmann even parameters a, b and Grassmann odd parameter c) 

{Tr(aA + bB + cx),Q e } = Tr(aSA + b5B + cSx) , (16) 

with 5A,5B,5x the supersymmetry variations given by Eq. (11) above, after elimination of 
the auxiliary fields F, G by their equations of motion. 

3. The Supersymmetric Yang-Mills Model 

As our next example of a component field supersymmetric model, we discuss super- 
symmetric Yang-Mills theory. (In Ref. [7] we have given a simpler analog of this discussion, 
in the context of the matrix model for M theory.) We start from the trace Lagrangian 

L = J £xTr[^F^ - xfD,* + \d 2 ] , (17a) 

with the field strength F^ v and covariant derivative constructed from the gauge potential 
A^ according to 

F^u = d/j,A u — d u A^ + [A^j A u ] 

D,0=d,0+[A^O] (176) 

^D^F uX + D„F x ^ + D x F^ = . 

In Eq. (17b), the potential components A^ are each an anti-self-adjoint, and the auxiliary 
field D a self-adjoint, bosonic N x N matrix (or operator), and each spinor component of 
X is a self- adjoint fermionic N x N matrix (or operator). The Euler-Lagrange equations of 
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motion are 

D =0 

7 "£> M X =0 (18a) 
D fl F^=2g 2 xYx ; 

as usual for a gauge system, the v — component of Eq. (18a) is not a dynamical evolution 
equation, but rather the constraint 

D e F eo = 2g 2 xi°X • (186) 
Going over to the Hamiltonian formalism, the canonical momenta are given by 

PA e = ~-^F oe , p x = ix T , (19a) 
and the axial gauge trace Hamiltonian (see [1] for a derivation and references) is 

H = H4 + H X , (196) 



with 



H.=/^Tr(^!^^-^F 2 3 



^ 2 {diA 2 - d 2 A x + [A^A 2 ]f - ^[(dsA.f + (d 3 A 2 f]) 



1 f°° 

F os =t;9 2 / dz'e(z - z') [-(p x x + X T P X ) + D 1 p Al + D 2 p M 

where we have taken care to write H in a form symmetric in the identical quantities p x and 
iX T , and where e(z) = 1(— 1) for z > 0(z < 0). The trace three momentum is 

P m = - / d 3 xTr F me p Ae +p x D m x) , (20a) 
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and the conserved operator C of Eq. (7) is given by 

C= [ d 3 x(J2[A e ,p Ae ]-{ x ,p x }) , (206) 
"* e=i 

with a contraction of the spinor indices in the final term of Eq. (20b) understood. By virtue 
of the constraint of Eq. (18b), the conserved operator C can also be written as 



2 

C = - I d 3 xJ2d ePAe = - [ d 2 S e p Ae , 

J J sphere at oo 



(20c) 



which vanishes when the surface integral in Eq. (20c) is zero. 
Making now the supersymmetry variations 

^=(^[7M,7,]^ + ^75£>)e (21a) 
5D ^ZisrfDpX , 

in the trace Lagrangian, we find that when e is constant, the variation vanishes. Again the 
calculation parallels that done in the c-number Lagrangian case, except that the trilinear 
cyclic identities of Eq. (lc) are used in place of commutativity/anticommutativity of the 
fields, and the vanishing of terms cubic in x is most easily established by using Eq. (12a) 
and the cyclic identity valid for Majorana representation 7 matrices (see the Appendix) 

E (?'fU(?'YvU = ■ (216) 

cycle a— >6— >d— >a 

When e is not a constant, the variation of L is given by 

5L = J d 3 xTr( J"9 M e) 



(22a) 



from which we construct the trace supercharge Q Q as 

Q Q = |d 3 xTr^(p x + ?X T )^^,7l« • (226) 
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Again, it is straightforward to check, using the equations of motion and the cyclic identity, 
that TrJ^ is a conserved trace supercurrent, which implies that the trace supercharge is 
conserved. 

We are now ready to check the closure of the supersymmetry algebra under the 
generalized Poisson bracket of Eq. (5a), which for the Hamiltonian dynamics of the super- 
symmetric Yang-Mills model gives 

^=*i±%W;-±%j£-(°"M ' (23a) 

We proceed now by the Fierz transformation method mentioned in Sec. 2 above. We begin 
by rewriting the boson terms in Eq. (23a) as 

2 

d 3 xTrJ2f>aP c g e c-{a^l3) , (236) 
and the fermion terms as 

r 4 

a <-> 0) , (23c) 

J d=i 

with the coefficient functions f,g,h,k readily determined once the operator variations in 
Eq. (23a) have been computed. Performing a Fierz transformation by using Eq. (A. 80) of 
West [6] then shows that verifying the supersymmetry algebra of Eq. (14b), with Q a ,p , H 
and P now given by the Yang-Mills expressions of this section, is equivalent to verifying the 
three identities 

(...v...) =iH 

(...y™...) =-P m (24a) 
2 



/■ 



(....[y.yi-) =0 

with 



(...r...) = - [ d 3 x\Tr (j2 f? T a9id + E hdar ak dc ) • (246) 
J \t=\ d=i J 
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Equation (24b) contains both local terms, and nonlocal terms that couple variables at dif- 
fering values of z. The nonlocal terms are found to vanish identically when the constraint 
equation and symmetries of the integrands are taken into account, while the local terms are 
seen, by an enumeration of cases, to obey Eq. (24a). Examining the role of the supercharge 
as a generator of transformations, in analogy with Eq. (16), in the Yang-Mills case the super- 
charge is found to generate the supersymmetry variations of Eq. (21a), plus an infinitesimal 
change of gauge. 

We conclude by showing how the results of this section include the conventional 
case of U(M) supersymmetric Yang-Mills theories, and why at the same time they are more 
general. Consider the case in which the N x N matrices acted on by Tr have dimension 
N = MP, and expand the matrices on a complete basis of U(M) matrices Aj, so that for 
the potential we have 

M 2 
i=0 

where the coefficients A^ are now P x P matrices. Then the commutator term LA^,^] in 
Eq. (17b) becomes 

ij 

When P — 1, so that the A^ all commute, the second term on the right hand side of 
Eq. (25b) vanishes, and it reduces to the conventional expression for the commutator term 
in a Yang-Mills theory. However, our formalism generalizes this conventional model to allow 
any P > 1, including the limit P — > oo, in which case the second term in Eq. (25b) contributes 
as well as the first. Similar remarks apply to the other matrix commutators appearing in 
the derivations of this section. 
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4. Superspace Considerations and Discussion 

The derivations of Sees. 2 and 3 have all been carried out in the component formal- 
ism, which requires doing a separate computation for each Poincare supersymmetry multi- 
plet. However, there is a simple and general superspace argument for the results we have 
obtained. Recall that superspace is constructed by introducing four fermionic coordinates 
9 a corresponding to the four space-time coordinates x^. The graded Poincare algebra is 
then represented by differential operators constructed from the superspace coordinates, and 
superfields are represented by finite polynomials in the fermionic coordinates 8 a , with coef- 
ficient functions that depend on x^. To generalize the superspace formulation to give trace 
dynamics models, one simply replaces these coefficient functions by N x N matrices (or oper- 
ators), and one inserts a trace Tr acting on the superspace integrals used to form the action. 
Then the standard argument that the action is invariant under superspace translations still 
holds for the trace action formed this way from the matrix components of the superfields. 
We immediately see from this argument why it is essential for the supersymmetry parame- 
ter e to be a c-number and not also a matrix; this parameter appears as the magnitude of 
an infinitesimal superspace translation, and since the superspace coordinates x M and 9 a are 
c-numbers, the parameter e must be one also. 

The construction just given gives reducible supersymmetry representations, and var- 
ious constraints must be applied to the superfields to pick out irreducible representations. 
Since these constraints act linearly on the expansion coefficients, they can all be immedi- 
ately generalized (with the usual replacement of complex conjugation for c-numbers by the 
adjoint) to the case in which the coefficient functions are matrices or operators. 

The simplicity of this argument suggests that for all nonextended rigid supersym- 
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metry theories for which there exists a superspace construction, there should exist trace 
dynamics generalizations, with component field forms analogous to those presented above 
and with corresponding representation covariant 7 matrix identities. In this paper we have 
not dealt with either extended supersymmetries, or with locally supersymmetric theories; 
these will be the subject of further investigations into supersymmetric trace dynamics theo- 
ries. 

Acknowledgments 

This work was supported in part by the Department of Energy under Grant #DE- 
FG02-90ER40542. I wish to thank Ed Witten, Andrew Millard, and members of the Prince- 
ton graduate student supersymmetry discussion group, for useful conversations. 

Appendix: Gamma Matrix Conventions and Identities 

We work with Majorana representation 7 matrices constructed explicitly as follows. 
Let 0X2,3 an d 71,2,3 be two independent sets of Pauli spin matrices; then we take 

7 = - 70 = -ia 2 Ti 
7° =H° = o 2 T\ 

7 1 =7i = °3 

7 2 =72 = -02T2 (A. la) 

7 3 =73 = -V\ 

75 =Z7 1 7 2 7 3 7° = -^7-3 

7° 75 =ir 2 , 
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so that 7°,75,7°75 are skew symmetric and 7 1 ,7 2 ,7 3 are symmetric, and 

f^ T f = _y . (A.lb) 

For this choice of 7 matrices, the four matrices 7^ are real. 

The identities of Eqs. (12b), (15), and (21b) are easily verified by the following 
method. Replace each four component spinor index a, b, c, d by a pair of two component 
spinor indices AA 1 ', BB', CO, DD', with the unprimed indices on the matrices a and the 
primed indices on the matrices r. Then each identity takes the form 

A-AA'BB'CC'DD' = , (A. 2a) 

with Eq. (A. 2) symmetric under the simultaneous interchange 

A <-> B, A' <-> B' , (A.2b) 

for the identities of Eq. (15) and (21b) (which are symmetric under a <-> 6) and antisymmetric 
under this interchange for the identity of Eq. (12b) (which is antisymmetric under a «-> b). 
To verify an identity, it suffices to verify the vanishing of its contraction with a complete set 
of 16 projectors on the 4x4 matrix with indices a, b, which in terms of a and r are 

[Sab, (01,2,3) ab] <8> V>a>b>, (^1,2,3) a< ,b>] , (-4-3) 

10 of which are symmetric and 6 of which are antisymmetric under the interchange of 
Eq. (A. 2b). Thus 10 contractions must be done to verify the identities of Eqs. (15) and 
(21b), and 6 contractions to verify the identity of Eq. (12b); these are readily done since 
the projectors involve two factors which repeat in different combinations, and since the 
contractions for individual factors involve only Pauli matrix arithmetic. 
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The identities of Eqs. (12b), (15), and (21b) are representation covariant, in that 
they do not take the same form in representations in which the Dirac gamma matrices are 
complex rather than real. To see this, we note that the matrices in a general representation 
7g are related to the Majorana representation matrices 7 M given above by 



with U a unitary matrix which in general is complex, as a result of which the row and column 
indices transform with different matrices. However, the identities used in the text mix row 
and column indices; for example, in Eq. (12b) there is one term in the cyclic sum in which 
o is a row index, and two terms in which a is a column index. (By way of contrast, the 
more familiar Fierz identities only interchange two row indices, and so do not mix row and 
column indices.) Hence we cannot get a representation invariant form of the identity by two 
applications of Eq. (A. 4), since in the second and third terms of the cyclic sum, we will have 
a row index contracted with a U and a column index contracted with a U*, which does not 
correspond with Eq. (A4). However, we can easily get a representation covariant form of 
Eq. (12b) by contracting all indices with a U*, and wherever U* contracts with a column 
index using the identity 



with 7 a matrix which appears on pp. 341-342 of the book cited in Ref. [1] (because it 
plays a role in the transformation properties of the Dirac equation in quaternionic quantum 
mechanics). We can then apply Eq. (A. 4) to all the gamma matrices, giving for Eq. (12b), 
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7 £ = ip^u = U T *^U 



(A4) 



u* = UU* T U* = UY , 7 = U T U 



(A5) 



for example, the representation covariant form 

E KtVWtV )cd + (7°757* U(7°757*) C J = • (A6) 

cycle a—*b^d^a 

For a change of representation which preserves reality of the 7 matrices, we have U* = 
U, 7 = U T U = U* T U = 1, and Eq. (A. 6) is identical to Eq. (12b), but for general changes 
of representation 

the identity is form covariant but not form invariant. 
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